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Recent results from systematic parameterizations of 
Ginsparg- Wilson fermions 
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^ Institut fiir Thcorctischc Physik, Universitat Regensburg, D-93040 Regensburg, Germany 

The "Fixed Point Dirac Operator" and "Chirally Improved Fermions" both use large numbers of gauge paths 
and the full Dirac structure to approximate a solution of the Ginsparg- Wilson equation. After a brief review of 
the two approaches we present recent results for quenched QCD with pion masses down to 210 MeV. We discuss 
the hmits and advantages of approximate parameterizations and outline future perspectives. 



1. Introduction 

The fundamental difficulties with the imple- 
mentation of chiral symmetry on the lattice 
were finahy overcome with the rediscovery of 
the Ginsparg- Wilson equation ||^,^ for the lattice 
Dirac operator D, 



D-f5 + i^D ^ 2aDR-izD 



(1) 



On the right-hand side i? is a local operator com- 
muting with 75 and a denotes the lattice spacing. 
This equation implies exact chiral symmetry on 
the lattice Results obtained with Dirac oper- 
ators obeying (|l|) exactly or approximately now 
allow to test QCD also in the chiral regime (see 
for a recent review) . The most widely used chiral 
Dirac operator is the so-called overlap operator 
|5|. The overlap formula gives a simple explicit 
prescription how to construct a Ginsparg- Wilson 
fermion (i.e. a solution of Eq. (|l])) from any de- 
cent lattice Dirac operator. An approach related 
to the overlap formalism are domain wall fermions 
|6) where an auxiliary fifth dimension is used to 
implement the chiral symmetry. 

Besides the overlap and domain wall fermions 
two more approaches to chiral symmetry are 
known, fixed point fermions [^j8P,|l^ and the 
chirally improved operator [ p[ 12| . The fixed 
point (FP) Dirac operator is constructed from 
the saddle point evaluation of the RG equations. 

*The work reported here was done in the BGR (Bcrn- 
Graz-Regensburg) collaboration. The author is supported 
by the Austrian Academy of Sciences, APART 654. 



The chirally improved (CI) operator is obtained 
by a systematic expansion of a solution of the 
Ginsparg- Wilson equation. In a practical applica- 
tion both the FP and the CI operator will use only 
finitely many terms (essentially restricted to the 
hypercube) and one can expect |ll3| only an ap- 
proximation of a Ginsparg- Wilson fermion. Ex- 
act chiral symmetry can, however, be obtained by 
using the FP or CI operators as a starting point 
in the overlap projection (see also ||l^). Using 
an improved operator for the overlap was found 
to improve the localization and one can hope to 
also obtain better dispersion and scaling proper- 
ties. However, additional overlap steps also drive 
up the cost of the Dirac operator in numerical 
implementations. 

In this contribution we analyze how far into 
the chiral region one can proceed with FP or CI 
fermions without additional overlap steps. Fol- 
lowing a short review of the construction of FP 
and CI fermions we will discuss several recently 
obtained results. We address the spectrum of the 
Dirac operator and the determination of the topo- 
logical charge and susceptibility from the index 
theorem. This is followed by a discussion of the 
results from quenched spectroscopy with empha- 
sis on the pion mass close to the chiral limit where 
we extract the quenched chiral log parameter 5. 
We study the scaling behavior of rho and proton 
masses and of the pion decay constant. We con- 
clude with a discussion of the range of pion masses 
where one can work competitively with FP and 
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Figure 1. Schematic representation of a general 
lattice Dirac operator. 

CI fermions without overlap projection. 

2. Construction of FP and CI fermions 

The construction for both the FP and CI op- 
erators starts from the expression for a general 
Dirac operator on the lattice |||,|ll|. In Fig. |l| we 
give a schematic representation of such a general 
lattice Dirac operator. 

Let us start the discussion of the structure with 
already familiar terms. In the dashed square 
in the upper left corner of Fig. |l| we show the 
terms that are used in the Wilson Dirac opera- 
tor. There are two terms coming with the unit 
matrix in Dirac space: Firstly a constant term 
consisting of the mass parameter and the cen- 
ter term in the discretization of the Laplacian. 
They are represented by a dot. Secondly also the 
nearest neighbor terms from the discretization of 
the Laplacian come with a unit matrix in Dirac 
space. These hopping terms are represented by 
single hops, i.e. straight lines in all directions (for 
artistic reasons we show only two of the four pos- 
sible directions). Besides the two terms trivial in 
Dirac space, the Wilson operator also contains 
the naively discretized massless Dirac operator 
where a sum over all 7^ together with hops in 
positive and negative /x-direction occurs (we show 
only one of the four directions). In order to de- 



scribe a derivative here the hops in positive and 
negative direction come with different signs. The 
parameters si, S2 and vi are real numbers chosen 
such that a fermion free of doublers with a given 
mass is described. 

A more general Dirac operator is obtained 
when more terms are used to discretize the deriva- 
tive, such as next to nearest neighbor terms, sta- 
ples etc. Similarly one can also extend the num- 
ber of terms in the trivial Dirac sector. Finally 
from the Symanzik improvement program it is al- 
ready known that also sectors of the Clifford alge- 
bra other than the trivial and vector sectors can 
contribute to the Dirac operator. So the struc- 
ture of a general Dirac operator as depicted in 
Fig. 1^ is a sum over all elements of the Clifford 
algebra. Each element is multiplied with paths 
of link variables and each path has some coef- 
ficient Si,Vi,ti,ai,pi. Some of the paths have 
the same coefficient but differ by relative sign 
factors. These sign factors are entirely deter- 
mined by implementing the symmetries C,P and 
75-hermiticity (i.e. 751575 = ZJt). Rotation in- 
variance requires paths related by rotations on 
the lattice to come with the same coefficient and 
translation invariance makes the coefficients in- 
dependent of the actual space time point. With 
our particular choice of the representation of the 
Clifford algebra the coefficients Si,Vi,ti, ai,pi are 
real. A further generalization can be obtained by 
allowing these coefficients to be real functions of 
local loops of gauge links. This option is used for 
the construction of the FP Dirac operator. 

The FP and the CI operator now differ only 
in the method for determining the coefficients 
Si,Vi,ti,ai,pi. For the FP operator the goal is 
to use the coefficients to approximate a solution 
of the fixed point equation for the Dirac operator. 
The basic idea is that near the critical surface the 
theory is invariant under real space renormaliza- 
tion group transformations that relate the quan- 
tum fields on a fine lattice to fields on a coarser 
lattice. In general the corresponding equations 
are quite involved but in the weak coupling limit 
they can be solved using the saddle point method. 
For the Dirac operator the FP equation reads 

Dc = K - K^n [Df + Kn^n]~^ n'' . (2) 
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Here Dc (Df) is the FP Dirac operator evalu- 
ated on the coarse (fine) gauge field configura- 
tions, n is the blocking kernel for the fermions 
and K a free parameter of the blocking proce- 
dure which can be used to maximize the slope 
of the exponential decay of the Dirac operator. 
The coefficients Si,Vi,ti,ai,pi of the Dirac oper- 
ator were now adjusted such that a functional 
measuring essentially the difference between the 
left-hand and right-hand sides of ^ on an en- 
semble of coarse and fine gauge configurations re- 
lated by the renormalization group is minimized. 
In the actual construction this procedure was it- 
erated starting from an ensemble at very weak 
coupling, followed by an intermediate step and 
finally a last minimization step at the target cou- 
pling. We remark that the action was optimized 
only for one such coupling corresponding to a lat- 
tice spacing of 0.16 fm and then used also for the 
other two couplings. The parameterized FP op- 
erator is described by 82 couplings corresponding 
to 41 independent terms and each of the coefR- 
cients Si,Vi,ti, Ui, pi is a linear function of local 
closed gauge loops. All terms of the FP operator 
are restricted to the hypercube. 

For the CI operator the strategy is to directly 
insert the expanded Dirac operator as depicted 
in Fig. 1^ into the Ginsparg- Wilson equation (|^) 
with a trivial kernel R = 1/2. The evaluation 
of the two sides of Eq. (|^) can be implemented in 
an algebraic computer program using directly the 
systematic representation of Fig. |^. The result 
is an expansion for both sides of the Ginsparg- 
Wilson equation similar to the expansion of the 
original Dirac operator. One can show that the 
individual terms (different element of the Clifford 
algebra, different paths) are linearly independent. 
The coefficients of equal terms on the two sides 
have to be equal and one can read off a system 
of coupled quadratic equations for the expansion 
coefRcients 

2si = sf + Ssj... + 8vl... 

2s2 = 2siS2 -I- 12S2S3... + 12wiW2... 

(3) 

This system is equivalent to the original 
Ginsparg- Wilson equation. After one truncates 



the expansion in Fig. |l| this system becomes fi- 
nite and can be solved numerically. In addition to 
the equations representing the Ginsparg- Wilson 
equation one can impose boundary conditions, 
i.e. additional equations for the coefficients en- 
forcing vanishing quark mass, the correct disper- 
sion relation in the free case, 0{a) improvement 
etc. The resulting CI operator is an approxima- 
tion of a solution of the Ginsparg- Wilson equa- 
tion. The CI operator has 19 coefficients and 
terms on the hypercube plus an extra L-shaped 
term of length \/5- 

3. Spectrum of eigenvalues 

A first impression of the quality of the approx- 
imation of a Ginsparg- Wilson Dirac operator can 
be obtained by inspecting eigenvalues of the Dirac 
operator in typical gauge backgrounds. For an ex- 
act solution of the Ginsparg- Wilson equation the 
spectrum is restricted to a circle of radius 1 with 
center 1 in the complex plane. For the two ap- 
proximations considered here the spectrum will 
not exactly lie on the Ginsparg- Wilson circle but 
show small fluctuations around it. 

In Fig. H we show a superposition of 6 spec- 
tra of the CI operator in quenched background 
gauge configurations on 12"^ x 24 lattices with 
lattice spacing a = 0.1 fm. Shown are the 50 
smallest eigenvalues for each configuration. It is 
obvious that the eigenvalues are not located ex- 
actly on the Ginsparg- Wilson circle but the fluc- 
tuations are rather small. In particular we do 
not find many configurations with negative real 
eigenvalues. Such so-called exceptional configura- 
tions give rise to numerical problems in the com- 
putation of the propagator and limit the small- 
est quark masses one can work at. As will be 
discussed below the suppression of the eigenvalue 
fluctuations achieved by the FP and CI operators 
allows us to work at considerably smaller quark 
masses than with e.g. Wilson fermions. 

An immediate consequence of the well ordered 
spectrum is the possibility to use the index the- 
orem ||l5| for evaluating the topological charge 
with the FP and CI operators, while for Wilson 
fermions the large fluctuations of the eigenvalues 
lead to a mixture of physical and doubler modes 
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Figure 2. Superposition of 6 spectra of the CI 
operator. Shown are the 50 smallest eigenvalues 
in the complex plane for 6 configurations on 12'^ x 
24 lattices with lattice spacing a = 0.1 fm. 



|lq| . We determined the topological charge as 
Q — — n+ where n_ ("•+) is the number of 
eigenvectors with negative (positive) matrix el- 
ement with 75. Measurements [p|jl7[| on differ- 
ent volumes V and lattice spacings give for the 
topological susceptibility {xtop — {Q'^)/y) val- 
ues of (196(4)MeV)'* for the FP operator and 
(191(2)MeV)4 for the CI operator. 

For the CI operator a detailed analysis of the 
behavior of Xtop across the finite temperature 
QCD phase transition was performed jT^. The 
quenched gauge configurations were generated 
with the Liischer-Weisz gauge action Q. Our 
results for the topological susceptibility as a func- 
tion of the temperature are displayed in Fig. |^. 
One finds that the results obtained on 6 x L'^ lat- 
tices with spatial extent L = 12, 16 and 20 nicely 
agree with each other showing that finite size ef- 
fects are under control. We also include results 
from 12^ and 16^ lattices to set the base line below 
Tc- The critical temperature Tc as determined for 
the Liischer-Weisz action in |l9) is marked by the 
dashed vertical line. One finds that the topologi- 
cal susceptibility starts to decrease already below 
Tc and quickly diminishes as the temperature in- 
creases further. The results obtained here using 
the index theorem agree well with calculations 
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Figure 3. The behavior of the topological sus- 
ceptibility across the QCD phase transition com- 
puted using the eigenvalues of the CI operator 
and the index theorem. The critical temperature 
is marked by a dashed vertical line. 



based on bosonic methods for the determination 
of the topological charge [|o| . 

4. Quenched spectroscopy 

During the last year the major goal of the BGR 
collaboration was to use the FP and CI Dirac op- 
erators for quenched spectroscopy. Before we dis- 
cuss the results of this study let us briefly outline 
the setting of these calculations. 

For the FP operator the perfect gauge action 
PH was used to generate the quenched ensembles. 
The gauge configurations were smoothened with 
perfect smearing ^ (which we consider as part of 
the parameterization of the Dirac operator) and 
subsequently fixed to Coulomb gauge. For the 
quark sources a Gaussian distribution was used. 

The CI operator was used with gauge ensem- 
bles generated with the Liischer-Weisz action . 
Here the gauge configurations were treated with 
one step of hypercubic blocking p2[. For the CI 
operator Jacobi smeared sources ]23| were used 
and no gauge fixing was necessary. 

For both operators we worked on three lattice 
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Figure 4. Pion mass (circles) and AWI mass 
(crosses) as a function of the bare quark mass. 
16^ X 32 lattice, a = 0.16 fm, FP operator. 



sizes, 8^ X 24, 12^ x 24 and 16^ x 32. For both 
gauge actions we used three different gauge cou- 
plings corresponding to lattice spacings of roughly 
a = 0.15 (0.16 for FP), 0.10 and 0.08 fm (a deter- 
mined from the Sommer parameter). The statis- 
tics varied between 100 and 200 configurations for 
each ensemble. Our choice for the volumes and 
lattice spacings allows to study the scaling be- 
havior with three different values of a at a fixed 
spatial length of about 1.3 fm and a finite size 
study at a fixed lattice spacing a — 0.16 fm. For a 
more detailed account of our setting for quenched 
spectroscopy see |Q . 

Let us begin the discussion of our quenched 
spectroscopy results with the pion mass and the 
axial Ward identity (AWI) mass. We computed 
the pion mass using different 2-point correlators 
(pseudoscalar, time component of the axial cur- 
rent, and the mixed correlator) and also com- 
pared correlators with smeared sink to point-like 
sink correlators. For all these measurements we 
find good agreement of the pion masses. The (un- 
renormalized) AWI mass ruAwi was computed as 

{doAoP) 



niAwi 



2(PP) 



(4) 



where P is the pseudoscalar density and Aq the 



time component of the axial vector current. In 
Fig. ^ we show the mass of the pion as computed 
from the pseudoscalar 2-point function (circles) 
and the AWI mass (crosses) as a function of the 
bare quark mass. This plot was generated using 
the FP operator on a 16^^ x 32 lattice at a lattice 
spacing a = 0.16 fm. 

It is a convincing sign of good chiral properties 
that both the pion mass and the AWI mass ex- 
trapolate to zero with the bare quark mass. How- 
ever, for smaller volumes ( L — 1.8 or 1.3 fm) we 
found that the quenched topological finite size ef- 
fects caused by the zero modes become more 
important. They can be removed by subtracting 
the scalar propagator, which has the same topo- 
logical finite size effect but larger mass, from the 
pseudoscalar 2-point function (see |24j for more 
details). In Fig. ^ we show also two fits to the pion 
data, a quadratic fit (full curve) and a fit includ- 
ing the quenched chiral log j26j (dashed curve). 
Our results for the quenched chiral log from a 
different method will be discussed below. 

Another test of chirality is to determine the 
residual quark mass. This was done using a lin- 
ear fit to the AWI mass and taking the intercept 
of the straight line with the horizontal axis as 
the residual quark mass rrires- For the FP op- 
erator we computed a residual mass of a rrires = 
—0.0006(4) at a = 0.16 fm increasing in size to 
arures = —0.0194(2) at a = 0.08. This increase 
is due to the fact that we used the FP action 
which was optimized for the a — 0.16 ensemble 
also at finer lattice spacing without redetermining 
the coefficients. For the CI operator the residual 
quark mass came out between arrires = 0.002(1) 
at a = 0.15 fm and arrires = 0.000(1) at a = 0.08 
fm. The smallest pion masses we have worked 
at are TOtt = 210 MeV for the FP operator and 
ttItt = 240 MeV for the CI operator. We expect 
that for both operators it is possible to go down 
to ~ 200 MeV without having to use exceed- 
ingly fine lattices. 

As another important benchmark measurement 
we looked at APE and Edinburgh plots. In Fig. |^ 
we show APE plots for the FP (top) and the CI 
(bottom) operator. For the FP operator we show 
three sets of data: The results for spatial extent 
L = 1.9 fm at a = 0.16 fm (circles), the same 
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Figure 6. Scaling of the octet baryon mass (cir- 
cles) and the vector mass (triangles). Filled sym- 
bols represent the FP results, while open symbols 
are used for the CI operator. 



Figure 5. APE plots for the FP (top) and CI 
(bottom) operators. 

ensemble but with the FP operator augmented 
with three steps of the overlap projection using 
Legendre polynomials (crosses) and finally results 
dX L — 2.6 fm, a — 0.16 fm (triangles). When 
comparing the FP operator on volumes with L = 
1.9 fm and L = 2.6 fm we find that the data 
agree well and we do not observe finite size effects 
for i > 1.9 fm. Also the overlap augmented FP 
operator (crosses) gives rise to results which are 
in very good agreement with the unimproved FP 
data. 

In the bottom figure for the CI operator we 
compare a data set from a finer but smaller lattice 
(a = 0.1 fm, L — 1.6 fm, circles) to results from 
a lattice similar to the one used for the FP oper- 
ator in the top figure (a = 0.15 fm, L = 2.4 fm, 
triangles). Here we do see a splitting between the 
two curves which we attribute to finite size effects 
and a small scaling violation. It is interesting to 
note that the data from the FP operator extrap- 
olate very well to the physical value (marked by a 
star in the plot), while the result from the CI op- 
erator undershoots the physical value similar to 



what is known from other quenched simulations 
(see e.g. 0). 

We conclude our discussion of the quenched 
spectroscopy with a brief discussion of the scaling 
properties of hadron masses. To do so, we work 
for both the FP and the CI operator at a bare 
quark mass which gives a ratio of mjr/wp = 0.7. 
At such a mass the statistical error is small and 
also a comparison with less chiral actions can be 
done (see e.g. j|]). In Fig. ^ we show the octet 
(circles) and vector (triangles) masses as a func- 
tion of the square of the lattice spacing. For both 
the horizontal and vertical axes we use the Som- 
mer parameter to set the scale. Filled symbols 
represent the FP results, while open symbols are 
used for the CI operator. The symbols are con- 
nected to guide the eye. 

The results from the two operators agree within 
error bars. Both sets of data show only a small 
deviation from a horizontal line indicating that 
both 0{a) and ©(a^) effects are small. Note that 
since both the parameterized FP operator as well 
as the CI operator are only approximate solutions 
of the Ginsparg- Wilson it can not a priori be ex- 
pected that 0{a) effects are absent. However, the 
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Figure 7. Pion decay constant as a function of 
the bare quark mass. 16"^ x 32 lattice, a = 0.1 fm, 
CI operator. 



plot shows that scaling violations are very small. 

5. Pion decay constant 

The definition of the pion decay constant can 
be combined with the axial Ward identity to yield 



2m \fZpp 



(5) 



Here Zpp is the prefactor of the pseudoscalar 2- 
point function at zero-momentum, 

Zpp 



{Pix,t)P{0,0)) 



2^m^ 



exp( — rriTri )• (6) 



Strictly speaking formula (|^) holds only for Dirac 
operators with exact chiral symmetry where 
the product of renormalization constants ZmZp 
equals 1. For the approximate solutions of the 
Ginsparg- Wilson equation which we use here this 
is currently only an assumption which will even- 
tually have to be tested. 

In Fig. 1^ we show according to Eq. as a 
function of the bare quark mass (both in units of 
the lattice spacing as determined from the Som- 
mer parameter) for the CI operator on the 16^ x 32 
ensemble at a = 0.1 fm. We interpolate the data 
using a second order polynomial which we subse- 
quently use to extrapolate the data to the chiral 
limit. We collected the results from the chirally 
extrapolated values of f^, from the CI operator 
and analyzed the scaling of this observable. We 



find that all data agree well within error bars and 
the discretization errors are small. Our quenched 
results come out 10-15 % larger than the experi- 
mental value (see for a more detailed discus- 
sion). 



6. Quenched chiral logs 

We have already briefly addressed the possibil- 
ity of extracting the quenched chiral log param- 
eter 5 from the pion mass (compare Fig. ^ and 
the discussion below the figure). However, here 
we use a method |^ which allows to increase 
the statistics by using pions with non-degenerate 
quark masses and which also gets rid of the depen- 
dence on the unknown chiral perturbation scale 

Ax- , , 

The idea is to form the combinations x and y 



mi +1712 ( m2 



■In , 

mi — m2 \mi 



y = 



4mim2 M- 



12 



(mi +7712)2 Mii^M22^ 



(7) 



Here M12 denotes the mass of the pseudoscalar 
with quark masses mi and m2. Using the results 
of [Eql one finds 



y 



1 + Sx + 0{m'^,5^,a^) , 



(8) 



i.e. to leading order the quenched chiral log pa- 
rameter can be read off from the slope in a plot 
of y versus x. 

In Fig. ^ we show such a x-y plot for the FP 
operator on the 16'^ x 32 lattice at a = 0.16 fm. 
The data lie inside a band which gives rise to 
S = 0.17(2) (compare (H] for previous determina- 
tions of S with chiral fermions) . For the CI opera- 
tor we obtain a very similar result of (5 = 0.18(3). 
These results were obtained by using the unrenor- 
malized AWI quark mass for mi,m2 in Eq. (0). 
For the CI operator we have also experimented 
with using the bare quark mass instead and we 
present a more detailed discussion in ||2^ . 

We remark that we also have obtained prelim- 
inary results for the pion scattering lengths and 
some details are also presented in [M. 



8 



0.95 - 
0.9 - 

0.85 - 
0.8 - 

0.75 - 
0.7 - 

0.65 — 



-1.5 



-0.5 



Figure 8. x-y plot used to determine the quenched 
chiral log parameter S (FP operator, 16"^ x 32, a ~ 
0.16 fm). 



7. Conclusions 

In this contribution we have presented re- 
sults from quenched QCD calculations using 
the parameterized FP and CI operators. Both 
these operators are approximate solutions of the 
Ginsparg- Wilson equation and are expected to 
have good chiral properties. Here these expec- 
tations are tested using various observables. 

We find that both the pion and the AWI mass 
nicely extrapolate to as a function of the bare 
quark mass. If one uses the offset in a linear ex- 
trapolation of the AWI mass as a criterion for 
the remaining chiral symmetry breaking we find 
residual quark masses between 1 and 4 MeV on 
a quite coarse lattice with a — 0.16 fm. We 
have successfully extracted the quenched chiral 
log parameter S using pions with non-degenerate 
quark masses. The octet and vector masses, as 
well as the pion decay constant show very good 
scaling behavior. We have demonstrated that we 
can reach pion masses of about 210 MeV with- 
out having to go to very small lattice spacings. A 
more detailed account of our measurements can 
be found in 
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Figure 9. A sketch comparing the cost of different 
Dirac operators as a function of the pion mass. 



We would like to wrap up our conclusions with 
a brief comment on where we think our approx- 
imate Ginsparg-Wilson fermions will be compet- 
itive. In Fig. y we show a sketch which com- 
pares the cost of different lattice Dirac operators 
as a function of the pion mass. For heavy quarks 
the operator of choice is certainly Wilson's Dirac 
operator, maybe combined with clover improve- 
ment. However, if one wants to go to smaller 
masses one has to increase the inverse gauge cou- 
pling to reduce eigenvalue fluctuations responsi- 
ble for exceptional configurations. This leads to 
a finer lattice and the volume has to be increased 
which drives up the cost, making Wilson fermions 
very expensive below pion masses of 300 MeV. 
At such small masses it is more economical to 
use the more expensive FP or CI operators which 
cost about a factor of 30 more than Wilson's op- 
erator (see ||2^ for a more detailed discussion of 
the cost). As already remarked, with these oper- 
ators we reached 210 MeV without having to go 
to very fine lattices and we expect that we will 
remain competitive down to at least 200 MeV. 
If one needs to go further into the chiral region, 
e.g. when computing the chiral condensate, one 
eventually has to use overlap projection to further 
decrease the pion mass. We believe that in the 
window between pion masses of 200 and 300 MeV 
the FP and CI operators are the best choice and 
we will explore whether this is sufficient to make 
reliable contact with chiral perturbation theory. 
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